Density of states near the Mott— Hubbard transition in the limit of large dimensions 
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The zero temperature Mott-Hubbard transition as a function of the Coulomb repulsion U is investigated 
in the limit of large dimensions. The behavior of the density of states near the transition at U = U c is an- 
alyzed in all orders of the skeleton expansion. It is shown that only two transition scenarios are consistent 
with the skeleton expansion for U <U C : (i) The Mott-Hubbard transition is "discontinuous" in the sense 
that in the density of states finite spectral weight is redistributed at U c . (ii) The transition occurs via a 
point at U = U c where the system is neither a Fermi liquid nor an insulator. 

PACS numbers: 71.30.+h, 71.27. +a, 71.28.+d 



The correlation-induced metal-insulator transition 
is a fundamental and challenging problem in condensed 
matter physics. Theoretical work has centered around 
the Hubbard model, where a metal-insulator transition is 
expected at half-filling for some critical on-site Coulomb 
repulsion U c . A new line of approach for studying the 
metal-insulator transition has recently been opened by 
the limit of high dimensions Q . At present this dynami- 
cal mean field theory seems the only tractable method to 
obtain exact statements concerning the Mott-Hubbard 
transition. Such results should be of considerable inter- 
est also for the physically relevant threedimensional case. 

In the dynamical mean field approach the behavior of 
the density of states pie) near the Mott-Hubbard transi- 
tion at half-filling and zero temperature is of particular 
interest. Georges et al. || have described in detail a tran- 
sition scenario where the spectral weight in the vicinity 
of the Fermi surface vanishes continuously as U t U c 
(compare Fig. 1): In some finite interval [ep — A, ep + A] 
around the Fermi surface one finds 
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This implies that at U c a finite excitation gap of size 2A 
opens in the density of states turning the Fermi liquid 
at once into an insulator. However, the analytical so- 
lution of the d — > oo Hubbard model is unfortunately 
still out of reach and the above transition scenario has 
been discussed controversially in the literature (see e.g. 
Refs. 0|). 

The purpose of this Letter is to establish constraints 
regarding possible transition scenarios, based on an anal- 
ysis of the dynamical mean field equations in all orders of 
the skeleton expansion. Conceptually the arguments are 
similar to the reasoning used by Luttinger || and Langer 
|jj in order to derive the (e — ep) 2 -behavior of the imag- 
inary part of the self-energy in the vicinity of the Fermi 
surface. An asymptotic sum rule for the imaginary part 
of the self-energy can be established that implies the fol- 
lowing constraint for the transition scenario: For all in- 



tervals with A > around the Fermi surface the limit 
in Eq. (Q) has to be nonzero as U | U c . Therefore the 
Mott-Hubbard transition investigated here is either (i) 
discontinuous in the sense that finite spectral weight is re- 
distributed at U e , or (ii) there is nonzero spectral weight 
in any neighborhood of the Fermi surface for U — U c . 

Possibility (ii) implies that the density of states is 
still gapless at the transition, i.e. pseudogap-like with 
p(ep) — but p(e) cx |e — £f| q , a > in the vicinity 
of ep. This would signal the existence of a non- Fermi liq- 
uid point separating metallic and insulating regimes. The 
analysis in this Letter will also provide insights in var- 
ious analytical approximation schemes devised to solve 
the dynamical mean field equations. 

The Hamiltonian of the Hubbard model is 
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where the hopping matrix elements are scaled as — > 
t/y/d to obtain a physically meaningful limit of large di- 
mensions 0J. In the sequel t sets the energy scale and 
half-filling corresponds to ep = 0. For d — > oo an ef- 
fective action describing the single-site dynamics of one 
fermionic degree of freedom (c 0q , Cq q ) can be written as 
(for details see Ref. ||) 
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All local correlation functions of the original Hubbard 
model can be derived from (g). The effective action 
is supplemented by a self-consistency condition relating 
the Weiss local field Q(t — t') and the local propagator 
G(t-t') = -{Tc (t)cI(t>)}. 

We will investigate the self-consistency problem by us- 
ing a Bethe lattice with large coordination number. On a 
Bethe lattice the self-consistency condition takes a par- 
ticularly simple form H 
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The discussion will be restricted to zero temperature, 
half-filling and the paramagnetic phase. Solutions with 
magnetic ordering are therefore excluded (^] . 

The fc-independent self-energy of the Hubbard model 
is given by E(e+) = - G- 1 (e+) or 
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where p(e) = — ilmG(e + ) is the local density of states 
and e + = e + i0 + . By introducing the principal value in- 
tegral A(e) = P duo one can split the self-energy 
into its real and imaginary parts E(e + ) = K(e) — iJ{e) 
leading to 



K(e) = e - t 2 A(e) - 
J(e) = -nt 2 p{e) + 
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Notice that p(e) is symmetric and therefore A(e) anti- 
symmetric at half-filling. 

We investigate solutions near the metal-insulator tran- 
sition coming from the metallic side U f U c M. A key 
role in this transition in the limit of large dimensions 
is played by the fixed value of the density of states at 
the Fermi surface p{cf) = l/nt for the half-filled case 
for U < U c jl(^] . This is due to the Fermi liquid property 
J(ep-) =0. In the zero temperature metal-insulator tran- 
sition scenario described by Georges et al. || the spectral 
weight w in the vicinity of the Fermi surface vanishes con- 
tinuously while still p(ep) — l/irt for U < U c . Rephrased 
mathematically, the local density of states separates in a 
low- and a high-energy part for U ] U c 



p(e) = pi(e)+p h (e) 
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where p/j(e) vanishes in some finite interval [—A, A] 
around the Fermi surface. The low-energy part pi(e) 
describes a quasiparticle resonance at the Fermi sur- 
face with a continuously vanishing spectral weight w as 
U t U c . Eq. (H) is identical to the ansatz of Moeller et 
al. [0, see also Fig. 1. 

Due to the pinning pi{ep) = l/irt the quasiparticle res- 
onance at the Fermi surface has a width of order wt. The 
existence of one single low-energy scale wt in this tran- 
sition scenario implies that the following scaling ansatz 
becomes possible in the limit w — * [O 
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The dimensionless function f(x) is normalized, 
dxf{x) = 1, and fulfills /(0) = 1/tt. This implies 



de p(e) = w -> for U T U c 
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FIG. 1. Hypothetical density of states near 
a metal-insulator transition with vanishing spectral weight 
in the vicinity of the Fermi surface. 



p/j(e) also has some w-dependence since it contains the 
remaining spectral weight 1 — w. This weak dependence 
will be of no importance in the following discussion. 

We now show that the above transition scenario cannot 
be realized within the skeleton expansion. It is taken 
for granted here that the transition occurs at a finite 
critical U c < oo as supported by numerical calculations 
for nonzero temperature (see Ref. ||). 

Before proceeding with actual calculations we present 
an intuitive argument why the scenario (0) is not con- 
sistent with the skeleton expansion. Consider a (hypo- 
thetical) density of states near a metal-insulator tran- 
sition as sketched in Fig. 1. For simplicity we assume 
that fix) has compact support in [—a, a]. Then for suf- 
ficiently small w the principal value integral A(e) devel- 
ops zeroes for e = ±eo with eo of order y/wt. Obviously 
awt <C eo < A and therefore eo lies in the gap p(±eo) = 
in the limit w — * 0. In this case A(eo) = leads to a 5- 
function-like contribution in J(e) according to Eqs. (Q) 
and ([|). Hence a density of states as depicted in Fig. 1 
enforces an imaginary part of the self-energy that van- 
ishes everywhere in the gap except for (^-functions at ±eo 
jl2j (notice that according to Eq. (^) J(e) is well-behaved 
for all other energies including the band edges). In the 
skeleton expansion the imaginary part of the self-energy 
is related to the available phase space for scattering pro- 
cesses. With the above ansatz for p(e), however, this 
phase space does not increase when passing from the 
small energy interval [—awt, awt] to [—A, A]. Therefore 
the resonance of J(e) in the gap cannot be explained in 
any order of the skeleton expansion. This argument can 
be generalized also to situations without a true gap. 

Let us formalize this reasoning. The method used is 
a demonstration by contradiction, i.e. we first assume 
that solutions of the d — > oo self-consistency conditions 
with the property (||) exist. In the sequel fix) need not 
have compact support. Consider the value of A(A) in the 
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limit w — > 0. The contribution from the states |e| < A in 
the principal value integral for A(A) behaves as 



A 



de 
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de p(e) = w/A , 



whereas the remaining spectral weight of order w for 
|e| > A contributes an essentially w-independent nega- 
tive term. Therefore eventually A(A) < 0, and thereby, 
according to (0) 



if (A) > A for w -> 



It will be demonstrated below that the skeleton expan- 
sion implies K(A) < 0, yielding a contradiction. In order 
to show this we derive a sum rule for the imaginary part 
of the self-energy. The self-energy can be expressed in 
the following manner in a skeleton expansion |l3| 

= [all possible skeleton diagrams with the 
unperturbed propagator Q(z) replaced 
by the true propagator G(z)] . (10) 

The diagrams in this series can be analyzed in an elegant 
way introduced by Langer 0: The imaginary part of 
a skeleton diagram is given by all possible Cutkosky- 
cuts across the internal lines. Diagrams with cuts across 
n internal hole and n+1 internal particle lines contribute 
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to J(e), where describes generalized vertex functions 
. In the vicinity of the Fermi surface this expansion 
gives the well-known behavior of the imaginary part of 
the self-energy for e — > ep M 
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where the coefficients T n are given by the above vertex 
functions and their derivatives in the infrared limit. The 
reason behind the increasing powers of e is the restriction 
of the available phase space for scattering processes in 
([□]) in the limit e -> 0. 

Next we investigate the behavior of f\ de J(e) in the 
limit w — > 0. The key observation is that the ansatz (||) 
automatically leads to a restriction of the available phase 
space for scattering processes on an energy scale smaller 
than A: We can consider the limit w — > in (|TT| ) and 
this leads to an expression for the integrated J(e) 
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The coefficients T n are averages of the functions 

r (n) (6> 6 - 6 - w 2 , • • ■ , w n +i) over tne Fermi liq- 

uid region: In ( |l l| ) the vertex functions are only probed 
on the Fermi liquid energy scale of order wt in the in- 
tegral J^ A de J(e). Therefore one can infer the scaling 
behavior of r„ with w from the behavior in the infrared 
limit: r„ oc T n for w — > |H. On the other hand from 



(9) the scaling ansatz (Q) one deduces r„aio 2n with an ' 



independent proportionality constant depending on the 
function f(x). Therefore the terms in the series ([l2]) are 
of order w leading to 
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de J(e) = awt 2 + 0(w 2 ) 
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with some dimensionless constant a that is finite for an 
integrable function f(x). Since K(e) and J(e) are con- 
nected by a Kramers-Kronig relation H 
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Eq. (O) implies that the positive contributions to K(A) 
in (p_4[) vanish for w — > 0, leading to 



if (A) < for w -> 



(15) 



Eqs. © and © are in obvious contradiction in the 
limit w — > 0: For sufficiently small but nonzero w, de- 
pending on the detailed structure of the upper and lower 
Hubbard bands and the function f(x), these two relations 
exclude each other. Therefore we obtain a contradiction 
already for some U < U c , that is in the metallic regime. 
This excludes solutions of the dynamical mean field equa- 
tions describing the transition scenario (Q) within the 
skeleton expansion. 

In this context a comment on the frequently used "iter- 
ated perturbation theory" (IPT) approach |L5) to d — > oo 
problems seems in order. IPT is second order pertur- 
bation theory for the self-energy in U where the Weiss 
propagator G(e + ) is used for the internal lines 

j( IPT >(e)=^/ 2 [ e dpp (e-p) 
Jo 

x / dvp (-v) p (p - v) (16) 
Jo 

with pa(e) — — ^ImQ(e + ). This approximation leads to a 
metal-insulator transition that is not seen when the full 
propagator is used for the internal lines in similar schemes 
Jl0| . Since the non-interacting density of states po{e) de- 
velops resonances on the energy scale ^/wt close to a hy- 
pothetical metal-insulator transition like in Fig. 1, sim- 
ilar resonances are found in j( IPT )(e). This mechanism 
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permits to fulfill the self-consistency conditions within 
the IPT-approximation. However, from the results in 
this Letter it is clear that such resonances in j( IPT )(e) 
are an artefact of the non-selfconsistent approximation 
used: If all orders of perturbation theory are summed 
up there is no phase space for scattering states available 
on the energy scale y/wt. This observation raises serious 
doubts whether IPT incorporates the correct mechanism 
that actually drives the transition. 

In summary in this Letter we have investigated the 
zero temperature metal-insulator transition in the half- 
filled Hubbard model in the paramagnetic phase on a 
Bethe lattice with large connectivity |Sj]. The framework 
used for this investigation was the skeleton expansion to 
all orders: The pointwise convergence of the skeleton ex- 
pansion in the metallic phase is the basic assumption (but 
compare fl6f| ) used in this Letter. Therefore our analysis 
was restricted to the metallic side of the transition since 
the skeleton expansion is known not to converge for an 
insulator. 

An important constraint for the transition scenario fol- 
lows from the asymptotic sum rule for the imaginary 
part of the self-energy. This sum rule leads to a competi- 
tion between the availability of phase space for scattering 
states and the suppression of spectral weight in the den- 
sity of states in the vicinity of the Fermi surface. This 
competition must be taken into account when investigat- 
ing the Mott-Hubbard transition in the limit of large di- 
mensions. It eliminates the scenario of a metal-insulator 
transition with vanishing spectral weight in the vicinity 
of the Fermi surface, compare Eq. (Q). Only two transi- 
tion scenarios are consistent with the skeleton expansion 
for U < U c : (i) A "discontinuous" transition occurs at 
U = U c in the sense that finite spectral weight is redis- 
tributed, (ii) For U — U c there is nonzero spectral weight 
in any neighborhood of the Fermi surface (pseudogap- 
behavior), corresponding to some non-Fermi liquid point 
separating metallic and insulating regimes. Careful nu- 
merical studies are required to establish whether such a 
pseudogap-solution of the self-consistency equations is 
possible. Finally, the competition effect demonstrated 
above can be expected to be of importance for nonzero 
temperature too. 

The author is greatly indebted to D. Vollhardt, 
P. G. J. van Dongen, F. Gebhard and W. Metzner for 
valuable discussions and many important remarks. Dis- 
cussions on related subjects with W. Hofstetter, J. Schlipf 
and M. Kollar are acknowledged. 
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